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ABSTRACT 

Multistage multirate signal processing is a qualifying technology that brings DSP techniques to the applications 
changing the rate of a signal. The Multirate filtering technique is widely used for meeting different sampling rates in the 
system. This paper envisages the implementation of multistage proposed filter structures of polyphase decomposition 
techniques with decimator and interpolator. This proposed filter structure is implemented using an analysis filter and 
synthesis filter with polyphase branches which reduces the computational complexity and adopts parallelism. The results 
analysis show in three different cases of proposed filter structure with multistage and compare with wavelet 9/7 filter to 
find the best removing noise elimination filter structure for ECG signal. 
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INTRODUCTION 

Wavelet transforms have been used in the field of signal and image processing. Recently, research on wavelet 
construction called lifting scheme, has been established by Wim Sweldens and Ingrid Daubechies [1]. 
This is also referred to as wavelet 9/7 filter. This scheme was developed on FIR-based discrete transform. In this, the input 
signal is fed into a low pass filter and high pass filter separately. The outputs of the two filters are then sub sampled. 
The original signal can be reconstructed by synthesis filters h and g which take the up-sampled Lowpass and high pass 
inputs. This scheme is also used in the Laurent polynomial representation of the filter and Euclidean algorithm. 
These schemes shows some limitations on sampling methods. 

Frequency converters sampling methods and narrow band filtering are known, allowing significant computational 
efficiency [2, 3]. However, current design procedures for these multistage and multilayer filters address the specification of 
each phase individually, rather than simultaneously optimizing all of the filter phases [4]. The authors [5] formulate an 
algorithm that optimizes multi-stage adaptive coefficients and also provides sufficient conditions for multi-track filter 
identification. A multilayer digital filter (MDF) is a digital filter that changes the input sampling frequency from the input 
signal to another signal. There are many applications in communication, image processing, digital audio and multimedia. 
In [6], the modern DSP system uses MDF with three factors. First of all, MDF is used in two digital systems with a 
different sampling rate. Second, MDF is the best approach to solving the complex filtering problem. Third, multilayer 
filtering is used in the construction of the multilayer filter bank. 
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Manually, constructing a multi-stage poly-phase filter (i.e., decimation (reduced sampling rate), 
filtered and polyphase interpolation (increased sampling frequency) is a time-consuming and high-risk process, which we 
focus on in this paper is the third factor.We developed the multistage multirate system support with multirate polyphase 
filters. This elimination process is easy for specific and computing a multistage multirate design. The main concept of 
proposed system, the input signal is decimated/interpolated and decomposed with M th polyphase filter branch with 
sampling rate conversion is called analysis filter in the first stage. In the second stage, the input signal is decomposed with 
Decimator and then the polyphase multistage filter by converting the sampling rate to the rational factor L/M that 
represents the test filter of the system proposed. This process is continued for another stage for getting more and more 
smooth signal. 


MATERIALS AND METHODS 

Database 

In this paper, we need to import ECG data for analysis, and so used the MIT-BIH Atrial Fibrillation database [7]. 
This database includes different types of noise such as power line noise, baseline noise, EMG noise, abrupt shift noise and 
electrosurgical noise. 


Proposed MultistageMultirate Polyphase Filter 


Figure 1 shows the general scheme of proposed multistage multirate system design. This system has been used the 
fractional sampling rate by cascading a factor-of-L interpolator with a factor-of-M decimator, where L and M are positive 
integer. We have used as H(z) decimation filter and G(z) as interpolation filter. These two filters operate with same or 
different sampling rates, they can be replaced with a single filter designed to avoid aliasing. The 1 st stage and 2 nd stage 
form an analysis filter and synthesis filter with polyphase branches. Configurations based on polyphase decomposition are 
convenient in applications where L and M are small numbers. Otherwise, filters of a very high order are requested. 
For proposed multistage Multirate filter system, design considers three cases which are described below. 


Input 

Signal 


Im 




1 st stage 



Nth stage 

Figure 1: Proposed diagram for Multistage Implementation of Sampling Rate Alteration System 


Case 1 


Consider the proposed design of low-pass FIR filter with poly-phase decomposition and decimation for ECG signal to 
meet the following specifications: 


Table 1 


Parameters 

Parameters Values 

Decimate factor (M) 

5 

Pass-band edge 

0.0871 

Stop-band edge 

2% 

Pass-band ripples 

40 dB 

Stop-bands ripples 

60 dB 


Impact Factor (JCC): 4.9564 


NAAS Rating 3.80 
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Designing the filter transfer function H(z) by Nth order Lowpass FIR digital filter from the above specifications 
and returns the filter coefficients in length as shown in Table I and its frequency response in Figure 2 where we 

take the sampling frequency. The transfer function is formed as 

H(Z) = h[0] 4- h[l]z -1 4 h[2] z -2 + h[3]z -3 4- h[4]z -4 + h[5] z -s 4- !i[6] z -6 4-h[7]z -7 -1- h[8] z -3 4 h[9]z -9 (1) 

Table 2: The Impulse Response h(n) or Filter Coefficients of Frequency Sampling Filter (N=9) 


n 

h[n] 

N 

h[n] 

0 

-0.0022 

9 

-0.0022 

1 

0.0669 

8 

0.0669 

2 

0.1097 

7 

0.1097 

3 

0.1610 

6 

0.1610 

4 

0.1928 

5 

0.1928 



Normalized Frequency (xn rad/sample) 


Figure 2: Frequency Response of Design Low Pass for FIR Filter 

Then performed decomposition of the filter transfer function into 5 poly-phase components using the 
equation (2). 


H(z} = EjU 1 ) -f z L E l (z s )+ z 3 E s (z s } + z a E 3 (z s ) + z 4 E 4 (z s ) 


( 2 ) 


Where, 

+- h[h] z" 5 

EjC*) = {fill] + h[S] x~ s ] 

F t (x) = {h[2] +h[7k- 5 ] 

E 3 (z) = {h[3] -Hhl^z- 1 ] 

E 4 (z) = {t[4] + Ji[9]z -S ] 

So, 

M—L 

H(jz) = T s~ k E k (z Mm ) 

(3) 

Where, 

[NiM: 

E<(z) = V MJfra + fclz-"0 <k < M- 1 

H=f? (4) 
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We set up decimate-by-M=5 for the input signal x[n]. In this step, the input {x[n]} is decomposed into the set of 5 

subsequences: {%2 [^]}> Here, for the cause {x [n]} we have x [-1] = 

x[-2] = x[-3] = x[-4] = 0. The 5 subsequences are filtered in the poly-phase branches and added together to give the 

decimated signal bfae l™]). The realization analysis filter structure of proposed system is shown in Figure 3. 



Figure 3: Proposed Multistage Multirate System for Case 1 


We chose interpolation-by- L=5 into the multistage multirate system structure using polyphase decomposion as 
shown in Figure 3. In this step, signal C^d^is used as an input to the interpreter.The signal {x[n]} = 1'^]} i s 

filtered in the parallel poly-phase branches and the set of 5 signals ,{<9o idi M }*{d2 (#4 W j? {#5 Ml, is 

obtained as the equations (8)-(12). The matrix G is composed of 5 row 
vectors: {#ol n ]}> {9il n \}> {92 M)’ {9 3 \9a M) The samples of the interpolated signal CVmtMlare stored column¬ 

wise in matrix U. The interpolated signal is considered as a synthesis filter is obtained simply by picking up the samples 
from matrix U. 


g, W = {h[0]r dffc (n> + b[5]y £at (n - 1}} ( 5 ) 

g L W = {h[l]ft«(n} + liE&ljWtn “ 1)} (6) 

Si W = {bfflXtar (n) + WiVtoc (n - 1}} (7) 

Si Cn> = {h[3] jfcrc (n) + Ji[S]> W (n - 1}} (8) 

St W = [b[4];y da , (n) + h[9] (n - 1>} ( 9 ) 


Taking inverse Z-transform of the equations (5) to (9), is obtained as 
={h[o]Y[ K ] + i![5lz- L M>]j 
G L tsr) = fh[l]Y[z] +Ji[6]z- l Y|>]3 
Gjtz) =pi[2]Y]z]-Hh[7]z- 1 Y[z]} 


Impact Factor (JCC): 4.9564 


NAAS Rating 3.80 
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G 3 Cz) = [h[3]Y[z] -f h[S]z -1 Y[z]} 


(13) 


G 4 Cz 3 =[h[4]YIz]4-h[9]z- 1 Y[z]} 


(14) 


Case 2 


In case 2, we considered a N th 


order low-pass digital FIR filter with poly-phase decomposition and 


decimation/interpolation for biomedical signals to meet the following specifications: 


Table 3 


Parameters 

Parameters Values 

Decimate factor (M) 

5 

Filter order, N 

81 

Cut-off frequency 

1 (M 


Table 4: The Impulse Response h(n) or Filter Coefficients of Frequency Sampling Filter (N=81) 


n 

h[n] 

n 

h[n] 

0 

0.00019 

81 

0.00019 

1 

-0.00020 

80 

-0.00020 

2 

-0.00057 

79 

-0.00057 

3 

-0.00079 

78 

-0.00079 

4 

-0.00072 

77 

-0.00072 

5 

-0.00032 

76 

-0.00032 

6 

0.00037 

75 

0.00037 

7 

0.00113 

74 

0.00113 

8 

0.00163 

73 

0.00163 

9 

0.00154 

72 

0.00154 

10 

0.00068 

71 

0.00068 

11 

-0.00079 

70 

-0.00079 

12 

-0.00240 

69 

-0.00240 

13 

-0.00342 

68 

-0.00342 

14 

-0.00317 

67 

-0.00317 

15 

-0.00138 

66 

-0.00138 

16 

0.00157 

65 

0.00157 

17 

0.00467 

64 

0.00467 

18 

0.00652 

63 

0.00652 

19 

0.00594 

62 

0.00594 

20 

0.00255 

61 

0.00255 

21 

-0.00286 

60 

-0.00286 

22 

-0.00840 

59 

-0.00840 

23 

-0.01161 

58 

-0.01161 

24 

-0.01050 

57 

-0.01050 

25 

-0.00448 

56 

-0.00448 

26 

0.00502 

55 

0.00502 

27 

0.01471 

54 

0.01471 

28 

0.02041 

53 

0.02041 

29 

0.01859 

52 

0.01859 

30 

0.00803 

51 

0.00803 

31 

-0.00913 

50 

-0.00913 

32 

-0.02742 

49 

-0.02742 

33 

-0.03929 

48 

-0.03929 
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Table 4 Contd., 

n 

h[n] 

n 

h[n] 

34 

-0.03741 

47 

-0.03741 

35 

-0.01717 

46 

-0.01717 

36 

0.02128 

45 

0.02128 

37 

0.07245 

44 

0.07245 

38 

0.12641 

43 

0.12641 

39 

0.17140 

42 

0.17140 

40 

0.19695 

41 

0.19695 


From the above sections, we found the transfer function H(z) by N th order Lowpass FIR digital filter and filter 
coefficients h[n] in length N+l as shown in Table 4. Its frequency response in Figure 4 where fs=500Hz. The transfer 
function is formed as 


H[z] = fr[0] 4 A[l]z 1 4 h[2]z 2 4- ■■■ 4 fr[81]z 81 (15) 



Figure 4: Frequency Response of Design Low Pass for FIR Filter 

Then performed decomposition of the filter transfer function into 2 poly-phase components using the equation 

(15). 

H(z) = E 1[f (z s } 4- *“ l E j > 3 ) (16) 

Where, 

E 0 = fr.[0] 4 h[2]z~ 2 4 h.[4\z~ 4 4- -■* -t- li[80]z _so 
Ei = ft[l]z _1 - 1 - fa[3]z -3 4fr[5]z -5 4 ■■■ 4 fit[81]z _B0 

E^q — fr[G] + h[6]z~ e 4 ft[12]z -12 -1 -b fr[78]z -78 

E 31 = fr[2]z -3 -F ft[ 8 ]z -8 4 h[14]z~ 14 4 -"4 h[80]z -80 

E 32 = h[4]z~ 4 4 fi[lQ]z -10 4 h[16]z~ Le H-h h[76]z _7& 

E 10 = tfllz " 1 + h\7]z~ 7 4 frri3lz" 13 4. . . - \-h\79]z~ 79 
E L1 = Ji[3]z -3 4 A[flz -9 4 Ml5]z -15 4. ... 4A[Sl]z”® L 
E l2 = fe[5]z -5s 4 A[ll]z” u 4 A[l7]z -L7 4 - 4 h[??]z- 77 

Decimated-by-M=2. In this step, the input M n ]iis broken down into the set of 2 subsequences: {^[m]}, 
[m]}. Here, for the causal {x[n]} we have x[-l] = x[-2] = 0. The 2 subsequences are interpolated by 3 in the poly-phase 
branches and added together to give the decimated signal ^da^l^^The actualization analysis filter structure of the 


Impact Factor (JCC): 4.9564 


NAAS Rating 3.80 
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proposed system is shown in Figure 5. 


For multistage purpose, further we have up-sampled & down-sampled the decimated signal 


m J 3 S everal 


times, such as, M=40 & L=40, M=20 & L=20, M=20 & L=40, M=40 & L=20. 


Case 3 



Figure 5: Proposed Multistage Multirate System for Case 2 


Consider the design specification of case 2, the decimated signal VVdM imjfof case 2 is used as an input to the 
interpolator and interpolated by L=2, then the 2 subsequences are decimated by 3 in the poly-phase branches and added 
together to give the interpolated signal [>W [m]}. The fulfillment analysis filter structure of proposed system is shown in 
Figure 6. 


The signal — tVs! at [™]}is filtered in the parallel poly-phase branches and the set of 2 signals (1-5': [n]}, 

fei, [n]} and the subsequences of this two signals {iJcu [n]}, {t&Ji [n]}, {feos [n]} { iSis- [n]}, j U?n [n]} {fei .2 [n]} is 


obtained as the equations (28)-(32) and added together to give the interpolated signal &W [m]}. 

{tfoM) = {g[0],g[M],g[2Ml 

kit”]} = {g[-l],g[M-l],g[2M-1],...} 

{&[«]} = {g[-2],g\M-2],g[2M-2],...} 

&«-U»]} = tel-M + 1] fl[l] fl[^+ n, a \2M + 1],...} (17) 

Where the equation can be represented based on the design specifications as, 

S,[n] = + ft[2]y 4}st [n - l] + M4]y dst [n - 2] + -■■ + A|B0]y dat [n - 40] (18) 

5 L [n] = ft[l]y d 4 t [n] +- -1] + h[5]y dst [ n - 2] + - + fclSl]y dat [ ,l - 40] ( 19 ) 

flw [«] = fc[0]y dat IrJ + tn - 3] + ft[l2]y d4f [n - 6] +- ■■■ -(-ft[73]y dsc [n - 39] ( 20 ) 

0 ,i.[n] = d2]y dat [n - 1] + MB]y dsf [n- 4] +- Ml4]y dsf [n - 7] +- - +- MB0]y da( [n -40] ( 21 ) 

S >2 [«] = k - 2] +■ Ml0]y dst [ti - 5] + fc[l0]y d3C [n — 3] +- ™ +- h[76]y dat [f. - 33] ( 22 ) 

flioM = fi[l]r ri ^[n] +A[7]'w,„ r [n - 3] + A[l3]T ri „dtt - - 39] (23) 

Sll [«] - ft™ - 1] + d9 [ti- 4] 4- ft[ 1S] y dsf [pi - "] 4- ■■■ -1- [ti - 40] (24) 
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8u [«] = MS]y dat h - 2] 4- A[ll] wk - 5] +A[l7]y^ f [n - S] + - 4- A[77]ft«[n - 3B 
Taking inverse Z-transform of the equations (18) to (25), is obtained as 
G. } lz] = Ji[0Mz] + hl2]z- l Y[z] + hl4]z~*Ylz] +■ - + M-30]z- 4O y[z] 

G L lz] = &[lMz] +■ i j[3]z- 1 Z!z] + h[S'.z~ 2 Ylz] f - + A [81] z^y lz] 

G w lz] = A[0]7[z] + A[6 ]z- 3 Z[z] + A[l2]z-*y[ir] + — + A[7S]z -w y[z] 
c fL [s] - ftl2]s- 1 rU] q- fe[a]i- 4 z[s] + Mi4]s-- zls] + - + 

Cj L [z]] = t[4]z 2 Y[z\ -t-ft[lO]z 5 r[z] -t-ft[ie]z °Ylz] + -t- h[75]z ™l r |>] 

G LO lz] = Ji[l]Z[z] + h[7]z~ 3 Y[z] -Hft[l3]z -6, Zfz] + + h[79]z~ 3 *Y[z] 

G L1 [z] = M3]z“ L yDr] + hl9]z~*Ylz] + MlSjz"- Ylz] +■ ■■■ +- A[Sl]r _4I? y[z] 

C L! [z] =A[5]z“ s y[z] +- h[ll]z~*Y\z\ + A[l7]z" # y[z] +■ - h\77]z~^Y[z\ 


(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


For multistage purpose, again the interpolated signal {tVtn t [m]} is decimate-by-M=2 and 2 subsequences are 
interpolated by 3 in the poly-phase branches as a synthesis filter and added together to give the signal JWt [m]. 



3[ -* 0 1 —hT) 
H G oi(^H T3M7 r7 l— 

— HSmfeM nHT 17 !—-t 

-3|—^ |z ‘1 —►( + > 



MU—&— 

—-RMtMT]—* 

Figure 6: Proposed Multistage Multirate System for Case 3 


RESULTS AND DISCUSSIONS 


To measure the performance of different noise removal methods, the distortion between original signal and 
reconstructed signal is measured by the mean square error (MSE), signal-to-noise ratio (SNR) and correlation coefficients 
values. But in this paper, the performance parameters such as SNR (in dB), MSE, correlation which will prove the 
robustness of these algorithms are not able to show in values because of lacking of pair to pair signal as shown in 
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Figure 7. Due to the process, proposed system have used different sampling rate for data conversion and filtering purposes. 

Figures 8 to 10 shows the result of a proposed multistage Multirate system using ECG signals in case 1-3 as 
considered in this paper. All case use multirate filter specification as shown in Table 2-4 on multistage multirate system for 
remove noisy ECG data. We have considered 1st stage as analysis and 2 nd stage synthesis filtering part is performed in 
Multirate system. For multistage operation, we considered repetition process. For case 1, proposed multistage multirate 
system is considered multirate poly-phase decomposition with decimation and interpolation. The sampling rate conversion 
factor is considered by down sampling factor M=5 and filter order N=9 for case I. Left side of column in Figure 8-10 (a) 
shows, output, smooth ECG signal with three different stages of proposed multistage Multirate system for case 1-3 and 
compared with wavelet 9/7 [1],[8] and their perspective frequency histogram. Figure 8 (a), shows better smooth signal 
where sampling rate is down by M=5 in first stage and the 2 nd stage increases the sampling rata by L=5 into the multirate 
polyphase filter. For better performance again, apply 3 rd stage as analyses filter and got smooth signal than 2 nd stage as 
shown in Figure 8 (a). For comparison of the proposed system and wavelet 9/7, we have used histogram analysis and 
autocorrelation function. Right side column of Figure 8 (a) shows the frequency distribution of noisy ECG signal, 
smoothed signal by proposed system and wavelet 9/7. The proposed systems have shown better frequency distribution 
shape where positions are slight because of sampling rate increase or decreases in the system as compared with wavelet 
9/7. 


4-4 — A——41—4'—4——Jl——4r— — 4--—4-- — 4 —4— 4— * 


Synthesis IBtcr )2»d xfge) 


1 -'J-/'— 

. ’ 

(e) Sample index 

Figure 7: Multistage Filtering of Noisy ECG Signal (MIT BIH AF 08215m.mat) for Case 1 


On the other sides in Figure 8 (b) shows a comparison based on autocorrelations. For Case 2, the sampling rate 
conversion factor is considered by down sampling factor M=2 and filter order N=81. Figure 9 (a), sampling rate is taken 
for rational sampling factor L/M=3/2 of ECG signals in the analysis filtering part, and show better signal than noisy signal. 
On the other side, the synthesis filter decreases the sampling rate by M= 10, and get a better signal than analyzed signal 
(1 st stage) as shown in the Figures 9 (a). For better performance again apply analyses filtering by L=10 on that synthesis 
filtered signal and got a smooth signal than 2 nd stage as shown in Figure 9 (a). For Case 3, in the analysis, filtering 
(1 st stage) part, sampling rate increased because of rational sampling factor L/M=3/2. The performance between the 
analysis filter and synthesis filter is also increases the sampling rate L/M=2/3 and gave better results than 1 st stage as 
shown in the Figure 10 (a). For better performance apply the output of the 2 nd stage into the synthesis filtered (3 rd stage) 
where the sampling rate was L/M=3/2. The results showed smooth signal than 2 nd stage. 


From the histogram of figure 8-10 (a) for case 1-3, it is also shown that stages 1-3 has a different frequency 
density with respect to their frequency distribution of various sampling rates and compared with wavelet 9/7. 
The histogram of 1st stage to 3 rd stage rates is bell-shaped with one peak for all cases. There are no gaps or extreme values. 
So, proposed systems give better results than wavelet 9/7 systems for every case. 
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Figure 8: (a) Multistage Filtering of ECG Signal & Frequency Response (b) Normalized Autocorrelation for Case 1 
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Figure 9: (a) Multistage Filtering of ECG Signal & Frequency Response (b) Normalized Autocorrelation for Case 2 

Now briefly justify histogram for only case 1, 1 st stage has a maximum value of 90 and a minimum value of -85 
with one peak between -35 and -60. 2 nd stage has a maximum value of 90 and a minimum value of -85, one peak between 
-35 and -60. 3 rd stage has a maximum value of 90 and a minimum value of -85 with one peak between -35 and -60. 
The wavelet 9/7 has skew left in the frequency distribution where a maximum value of 140 and a minimum value of -110 
with one peak between 18 and -12. Such as case 2 & case 3 also show better smooth ECG signal for 1 st to 3 rd stage than 
wavelet 9/7. 


* -— v|_^— a| 

5 ~‘kl — 


%4~lr^r 


'• ‘ [a-I- -~v^----u|—--*j 

(a) 



(b) 


Figure 10(a): Multistage Filtering of Noisy ECG Signal (MIT BIH AF 08215m.mat) & Frequency Distribution by 

Histogram (b) Normalized Autocorrelation for Case 3 

For comparison between the ECG signal and proposed system, we have also considered autocorrelation processes 
for robustness of algorithms. Figure 8-10(b), the representation of auto correlation values of ECG signal as red color, 
1 st , 2 nd , 3 rd stage as green, blue black color and wavelet 9/7 as cyan color. The results of cases 1-3 with three different 
stages show high peaks and good correlated values with different sampling rate. Figure 8(b) for case 1 it also observed that 
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the autocorrelation of 1 st and 3 rd stage overlap to each other because the output of 1 st stage is applied to the 2 nd stage. 
Then the output of the 2nd stage is applied to the 3 rd stage and finally the L/M is same for l st & 3 rd stage that’s why their 
frequency distribution also same. 

CONCLUSIONS 

The novel multistage Multirate system has been proposed for removing noise,requiring low-cost and equal or 
fractional sampling rates for ECG signal processing. This system has considered three cases for better justification of the 
proposed system with sampling rate conversion. From the results analysis, it can be concluded that all cases show better 
quality of ECG signal with the effect of reduce/increase the sample values of noisy ECG signals. These algorithms can be 
more efficient for diagnosis purpose. In future, we will develop more multistage methodology for biomedical applications 
to get more accurate performance. 
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